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The statistical mechanics of flexible surfaces with internal elasticity and shape fluctuations is 
summarized. Phantom and self-avoiding isotropic and anisotropic membranes are discussed, with 
emphasis on the universal negative Poisson ratio common to the low-temperature phase of phantom 
membranes and all strictly self-avoiding membranes in the absence of attractive interactions. The 
study of crystalline order on the frozen surface of spherical membranes is also treated. 

I. INTRODUCTION 

The statistical mechanics of polymers, which are one- dimensional chains or loops to a first approximation, has 
proven to be a rich and fascinating field. P, Q, ll| The success of physical methods applied to polymers relies on 
universality - many of the macroscopic length scale properties of polymers are independent of microscopic details 
such as the chemical identity of the monomers and their interaction potential, j^j 

Membranes are two-dimensional {2D) generalizations of polymers. The generalization of polymer statistical me- 
chanics to membranes, surfaces fluctuating in three dimensions, has proven to be very rich because of the richer 
spectrum of shape and elastic deformations available. In contrast to polymers, there are distinct universality classes 
of membranes distinguished by the nature of their short-range order. There are crystalline, fluid and hexatic membrane 
analogues of the corresponding phases of strictly two-dimensional systems (monolayers) where shape fluctuations are 
frozen. 11110 




FIG. 1: Esfera (Sphere) 1976: Gertrude Goldschmidt (Gego). Stainless steel wire - 97 x 88 cm (Patricia Phelps de Cisneros 
Collection, Caracas, Venezuela). 

The closest membrane analogue to a polymer is a 2D flshnet-like mesh of nodes with a fixed coordination number 
for each node. A fixed-connectivity membrane with spherical topology from the world of art||^ is shown in Fig^ 
Bonds are assumed to be unbreakable while the nodes themselves live in flat d-dimensional Euclidean space R'', with 
a physical membrane corresponding to the case d = 3. The intrinsic crystalline order of fixed-connectivity membranes 
with, say, typical coordination number 6, leads to the alternative terminology crystalline membranes. They are also 
referred to as polymerized or tethered membranes. In general the Hamiltonian for a fixed-connectivity membrane 
will include both intrinsic elastic contributions (compression and shear) and shape (bending) contributions, since the 
membrane undergoes both types of deformation. _10J 

Flexible membranes are an important member of the enormous class of soft condensed matter systems, [sl ITU fl^ IT^ 
Soft matter responds easily to external forces and has physical properties that are often dominated by the entropy of 
thermal or other statistical fluctuations. 

This chapter will describe the properties of flxed-connectivity membranes with focus on the universal negative 
Poisson ratio that illustrates the novel elastic behavior of the extended (flat) phase of physical membranes, the 
tubular phase of anisotropic membranes and ordering on frozen curved membrane topographies. 
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II. PHYSICAL EXAMPLES OF MEMBRANES 



One can polymerize suitable chiral oligomeric precursors to form molecular sheets. [T^ This approach is based 
directly on the idea of creating an intrinsically two-dimensional polymer. Alternatively one can permanently cross- 
link fluid-like Langmuir-Blodeett films or amphiphilic layers by adding certain functional groups to the hydrocarbon 



dget t ni: 

tails and/or the polar heads, as shown schematically in Fig[21 

Lujuid membrane .... 



Polymerized membrane 

FIG. 2: Making a fixed-connectivity membrane by polymerizing a fluid membrane. 




FIG. 3: An electron micrograph of a 0.5 micron square region of a red blood cell cytoskeleton at magnification 365,000:1. The 
skeleton is negatively stained and has been artificially spread to a surface area nine to ten times as great as in the native 
membrane. Image courtesy of Daniel Branton (Dept. of Biology, Harvard University). 

The 2D-cytoskeletons of certain cell membranes are beautiful and naturally occurring fixed-connectivity membranes 
that are essential to the function and stability of the cell as a whole. ^ 17, .18. .19- ..2Q. JlJ The simplest and most 
thoroughly studied example is the cytoskeleton of mammalian erythrocytes (red blood cells). The human body has 
roughly 5 x 10^^ red blood cells. The red blood cell cytoskeleton is a protein network whose links are spectrin 
tetramers (of length approximately 200 nm) meeting at junctions composed of short actin filaments (of length 37 
nm and typically 13 actin monomers long)|2l|2l|23 (see FigOland Figg. There are roughly 70,000 triangular 
faces in the entire mesh which is bound as a whole by ankyrin and other proteins to the cytoplasmic side of the 
other key component of the cell membrane, the fluid phospholipid bilayer. Without the cytoskeleton the lipid bilayer 
would disintegrate into a thousand little vesicles and certainly the red blood cell would not be capable of the shape 
deformations required to squeeze through narrow capillaries. 

There are also inorganic realizations of fixed-connectivity membranes. Graphitic oxide (GO) membranes are micron 
size sheets of solid carbon, with thicknesses on the order of lOA, formed by exfoliating carbon with a strong oxidizing 
agent. Their structure in an aqueous suspension has been examined by several groups. ,25, 26, 27] Metal dichalcogenides 
such as M0S2 have also been observed to form rag- like sheets. j28j Finally similar structures occur in the large sheet 
molecules, shown in FiglSj and believed to be an ingredient in glassy B2O3. 

III. PHASE DIAGRAMS 

Let us consider the general class of _D-dimensional elastic and fiexible manifolds fluctuating in d-dimensional Eu- 
clidean space. Such manifolds are described by a d-dimensional vector r(x), where x labels the D-dimensional internal 



FIG. 4: An extended view of the spectrin/actin network which forms the cytoskeleton of the red blood cell membrane. Image 
courtesy of Daniel Branton. 




FIG. 5: The sheet molecule 8203. 



coordinates, as illustrated in FigEl A physical membrane, of course, corresponds to the case {D = 2,d = 3). 




FIG. 6: The parametrization of a membrane with internal coordinates x and bulk coordinates r(x). 



The Landau free energy of a membrane must be invariant under global translations, so the order parameter is given 
by derivatives of the embedding r, viz. the tangent vectors ta — -§^, with a = 1. • • • , D. Invariance under rotations 
in both the internal and bulk space limits the Landau free energy to the form 12^ l30| 



F{r) 



+ - / d^xd^y<5''(f(x)-f(y)) , 



(1) 



where higher order terms are irrelevant in the long wavelength limit. The physics of Eq.JQ depends on the elastic 
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moduli t, u and v, the bending rigidity k and the strength of self-avoidance b. The limit 6 = describes a phantom 
membrane that may self-intersect with no energy cost. 

For small deformations from a reference ground state one may write r(x) as 



f(x) = (Cx + u(x),Mx)) , (2) 

where u(x) are D "internal" phonon modes and /i(x) d — D "out-of-plane" height fluctuations. The case ( — 
corresponds to a mean field isotropic crumpled phase for which typical equilibrium membrane configurations have 
fractal Hausdorff dimension dn {dn = oo for phantom membranes) and there is no distinction between the internal 
phonons and the height modes. The crumpled phase is illustrated in FigCfa). 

The regime C 7^ describes a membrane which is "flat" up to small fluctuations. The full rotational symmetry of 
the free energy is spontaneously broken. The fields h are the Goldstone modes and scale differently than the phonon 
fields u. Figl^c) is a visualization of a typical configuration in the "flat" phase. 
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FIG. 7: Typical configurations of phantom membranes: (a) the crumpled phase, (b) the critical crumpling phase and (c) the 
flat or bulk-orientationally-ordered phase. Images are from the simulations of Ref.31. 

Phantom membranes are by far the easiest to treat analytically and numerically. They may even be physically 
realizable by synthesizing membranes from strands that cut and repair themselves on a sufficiently short time scale 
that they access self-intersecting configurations. One can also view the analysis of the phantom membrane as the first 
step in understanding the more physical self-avoiding membrane. Combined analytical and numerical studies have 
yielded a thorough understanding of the phase diagram of phantom fixed-connectivity membranes. 



A. Phantom Membranes 



The phantom membrane free energy is 



F{r) = J d^x 

The mean field effective potential, using the expansion of Eq.(j21), is 



(3) 



(4) 
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with minima 





t 

' 4:{u+vD) 



t > 
t < 



(5) 



This imphes a "flat" (extended) phase for t < and a crumpled phase for t > 0, separated by a continuous crumpling 
transition at t = 0, as sketched in FiglHl 




t< 




flat phase 




/ >0 




crumpled phase 



FIG. 8: The mean field free energy density / of fixed-connectivity membranes as a function of the order parameter t together 
with a schematic of the low temperature flat ordered phase and the high temperature crumpled disordered phase. 

Of course anything is possible in mean field theory but a variety of analytic and numerical calculations indicates the 
true phase diagram of phantom membranes is qualitatively like Fig|51 9] The crumpled phase is described by a line of 
equivalent Gaussian fixed points (GFPs). There is a crumpling transition line in the n — t plane containing an infrared 
stable fixed point (CTFP) which describes the long wavelength properties of the crumpling transition. Finally, for 
large enough values of k and negative values of t, the system is in a "flat" phase whose properties are dictated by an 
infrared stable flat phase fixed point (FLFP). 



crumpling transition line 



1/K 

crumpled phase 




FIG. 9: Schematic plot of the phase diagram for phantom membranes. 

1. The crumpled phase 
In the crumpled phase, the free energy Eq.©, for D >2, simplifies to 



F{r) = — J d^x {dav)'^ + irrelevant terms , 



(6) 
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since the model is completely equivalent to a linear sigma model with 0{d) internal symmetry in D > 2 dimensions 
and therefore all derivative operators in f are irrelevant by power counting. [32i. j33j The parameter t labels equivalent 
gaussian fixed points, as depicted in FigEl In renormalization group language there is a marginal direction for positive 
t. The large distance properties of this phase are described by simple gaussian fixed points with the exact connected 
2-point function: 



{ log |x| -.0 = 2. 

The associated critical exponents may also be computed exactly. The Hausdorff (fractal) dimension dn, or cquivalcntly 
the size exponent v ~ D/du, is given (for the physical case D — 2)hy 



= oo (i^ = 0) ^ i?G - logL , (8) 

where i?g is the radius of gyration and L is the linear membrane size. This result is confirmed by numerical simulations 
of fixed-connectivit y rn cmbranes in the crumpled phase where the logarithmic behavior of the radius of gyration is 
accurately checked. |3H-(34j-p5l| 



2. The Crumpling Transition 



Near the crumpling transition the membrane free energy is given by 



F{r) = j d^x 



(9) 



where the bending rigidity has been scaled out and v = v — j^. By naive power counting the directions defined by the 
couplings u and v are relevant for D < A and the model is amenable to an e = 4 — D-expansion. The /3 functions are 
given by0 



Pu{uR, vr) = -euR + ^ I + Y^j "fl + ^URVR + '^v\ \ (10) 



and 



I3v{ur,vr) = -£VR+ -^iy—u\ + —URVR + {Ad + b)v%^ . (11) 

These two coupled beta functions have a fixed point only for d > 219. This suggests that the crumpling transition 
is first order for d = 3. Other analyses, however, indicate a continuous crumpling transition. A revealing extreme 
limit of membranes was studied by David and Guitter. |46j| This is the limit of infinite elastic constants in the flat 
phase. Since the elastic terms in the Hamiltonian scale like in momentum space, as compared to q'^ for the bending 
energy, this limit exposes the dominant infrared behavior of the membrane. In this "stretchless" limit the elastic 
strain tensor Uap must vanish and the Hamiltonian is constrained, very much in analogy to a nonlinear sigma model. 
The Hamiltonian becomes 



Hnl^ J d^'a^iAr)^ , (12) 



together with the constraint dafdfjr = Sap- Remarkably, the /3-function for the inverse bending rigidity a — l/n may 
be computed within a large-d expansion, giving 

^, , da 2 /I const. \ , , 
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For d = oo there is no stable fixed point and the membrane is always crumpled. To next order in 1/d, however, 
Ea. (|13|l reveals an ultraviolet stable fixed point at a = Sir/d, corresponding to a continuous crumpling transition. 
The size exponent at the transition is found to be^47j 



2d 



1 



1 - - = 2/3 (for d = 3) 



(14) 



Le Doussal and Radzihovsky^^ analyzed the Schwinger-Dyson equations for the model of Eq.0 keeping up to four 
point vertices. The result for the Hausdorff dimension and size exponent is 



rfff = 2.73 



v = 0.73 . 



(15) 



Finally Monte Carlo Renormalization Group simulations |49l| of the crumpling transition find a continuous transition 
with exponents 



c\ ryi-! f 1 r-\ \ 



n m /n\ 



(16) 



Thus three independent 



iize exponent. 




(l + K) 



FIG. 10: Plot of the specific heat observable from the simulations of Ref.31. The growth of the specific heat peak with system 
size indicates a continuous transition. 

Further evidence for the crumpling transition being continuous is provided by numerous numerical simulations 

lisf where the analysis of observables like the specific heat (see Fig llO|l or the radius of gyration radius give textbook 

continuous phase transitions, although the value of the exponents at the transition are difficult to determine precisely. 



3. The Flat Phase 




FIG. 11: Membrane coordinates appropriate for analyzing fluctuations in the flat phase. 
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In a flat membrane (see Fig lll|l . it is natural to introduce the strain tensor 



Mq/3 = daUp + dpUa + dahdph 
The free energy Eq.ljSJ in these variables becomes 



(17) 



F(u, h) 



(18) 



where irrelevant hig her derivative terms have been dropped. One recognizes the standard Landau free energy of 
elasticity theorv. |5(]l| with Lame coefficients /i and A, plus an extrinsic curvature term, with bending rigidity k. These 
couplings are related to the original ones in Eq.© by /i = A = 2vC,'^~^ , k = kC,^~^ and t — — 4(/i+ -jA)C^~^- 

The large distance properties of the flat phase for flxed-connectivity membranes are completely described by the free 
energy of Eq. (|18|l . Since the bending rigidity may be scaled out at the crumpling transition, the free energy becomes 
a function of and The /3- functions for the couplings /i and A in the e-expansion are |5ll Is^ 



Afl) 



-sXr + ^ I y^l + 2(4 + 10A)Xrhr + 2dcXji ] , 



(19) 



1 /dc 



where is the codimension d — D and A = 2jf^p^Ti^- These coupled /3- functions possess four flxed points (see Fig ll2|l 
whose values are shown in Table HI 




FIG. 12: Phase diagram for the phantom flat phase. There are three infra-red unstable fixed points, labelled by FLPl, FLP2 
and FLP3, but the physics of the flat phase is governed by the infra-red stable fixed point (FLFP). 

The phase diagram revealed by the e-expansion is thus a little more complex than that sketched in Fig|5| The three 
additional flxed points are infra-red unstable, however, and can only be reached for very specific values of the Lame 
coefficients. 



4- The properties of the flat phase 



FigCJc) shows a typical equilibrium configuration for a membrane that has developed a preferred orientation in 
the bulk - the surface normals clearly have long-range order. In this phase the membrane is a rough extended 
two-dimensional structure. The rotational symmetry of the full free energy is spontaneously broken from 0{d) to 
0{d — D) X 0{D). The remnant rotational symmetry is realized in Ea. H18|l as 
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/i,(x) ^ /i,(x)+A*"x„ ; (20) 

where A^" is a D x {d ~ D) matrix. This relation provides Ward identities which greatly simplify the renormalization 
of the theory. 

The phonon and height propagators in the infrared limit are given by 

r™(p) ^ bf ; (21) 

where the last equation defines the anomalous bending rigidity k(p) ~ |p| The two scaling exponents rju and t] are 
related by the scaling relation 'Fl^ 

?7„ = 4 - - 277 , (22) 

which follows from the Ward identities fEa. H2U|) ') associated with the remnant rotational symmetry . The roughness 
exponent ^, which measures the growth with system size of the rms height fluctuations transverse to the flat directions, 
is determined from 77 by the further scaling relation 

C = ^ ■ (23) 

The long wavelength properties of the flat phase are described by the FLFP (see Figll2|). Since the FLFP occurs 
at non-zero renormalized values of the Lame coefficients, the associated critical exponents are clearly non-Gaussian. 
These key critical exponents have also been determined by independent methods. 

Large scale simulation of membranes in the flat phase model were performed in Ref.[23| The results obtained for the 
critical exponents are very accurate: 



77„ = 0.50(1) ; 77 = 0.750(5) ; C = 0.64(2) . (24) 



A review of numerical results may be found in R,efs.|44L l45|. 

The SCSA approximation gives a beautiful result for general d: 

V{d) = , ^ , (25) 

4 + \/16 - 24 + dl 

which for d = 3 gives 

77„ = 0.358 ; 77 = 0.821 ; C = 0.59 . (26) 

Finally the large-d expansion^^ gives 

1=\=^ 77(3) = 2/3 . (27) 



TABLE I: The fixed points and critical exponents of the flat phase. 
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The numerical simulations are in qualitative agreement with both the SCSA and large-d analytical estimates. 
On the experimental side we are fortunate to have two measurements of the key critical exponents for the flat phase 
of fixed-connectivity membranes. The static structure factor of the red blood cell cytoskeleton has been measured 
by small-angle x-ray and light scattering, yielding a roughness exponent of C = 0.65(10). '24j Freeze-fracture electron 
microscopy and static light scattering of the conformations of graphitic oxide sheets reveal flat sheets with a fractal 
dimension cIh = 2.15(6). Both these measured values are in good agreement with the best analytic and numerical 
predictions, but the errors are still too large to discriminate between different analytic calculations and to accurately 
substantiate the numerical simulations. 

The Poisson ratio |50|| of a phantom fixed-connectivity membrane (which measures the transverse elongation due to a 
longitudinal stress) is universal and within the SCSA approximation is given by 



<^iD) = =^ a(2) = -1/3 . (28) 

This result has also been checked in numerical simulations, [sj |3 Rather remarkably, it turns out to be negative. 
While Ref.jssjl finds a « —0.15 the latter simulation finds a ~ —0.32. Materials with a negative Poisson ratio 
have been dubbed auxetics\5^ . The wide variety of potential applications of auxetic materials suggests a fascinating 
role for flexible fixed-connectivity membranes in materials science (see Sect. Hvj). 

A final critical regime of a flat membrane is achieved by subjecting the membrane to external tension. |4 71] This gives 
rise to a low temperature phase in which the membrane has a domain structure, with distinct domains corresponding 
to flat phases with different bulk orientations. This describes, physically, a buckled membrane whose equilibrium 
shape is no longer planar. 



B. Self-avoiding Membranes 



Physically realistic fixed-connectivity membranes will have large energy barriers to self-intersection. That is they 
will generally be self-avoiding. Self-avoidance is familiar in the physics of polymer chains and may be treated by 
including the Edwards- type delta- function repulsion of the Hamiltonian in Eq.(^. A detailed summary of our current 
understanding is given in Kefs.!, [13 The essential finding is that self-avoidance eliminates all but the flat phase. 



1. Numerical simulations 



Numerical simulations are currently essential in understanding the statistical mechanics of self-avoiding membranes 
because the treatment of nonlinear elasticity together with non-local self-avoidance is currently beyond the realm of 
analytic techniques. 

Two discretizations of membranes have been adopted to incorporate self-avoidance. The balls and springs class of 
models begins with a network of N particles in a intrinsically triangular array and interacting via a nearest-neighbor 
elastic potential 



V..ir^ = \ ' f I'^f^ , (29) 
I oo tor |r| > 

where the free parameter b plays the role of a tethering length. An additional hard sphere steric repulsion forbids any 
node to be closer than a distance a from any other node: 



I lor |r| > (7 

Early simulations [3 [s^ of this class of model gave a first estimate of the fractal dimension for physical membranes 
compatible with the Flory estimate dn = 2{d + D)/{2 + D) = 2.5. [s^ The system sizes simulated, however, were quite 
small and subsequent simulations for larger systems found that the membrane is flat. |56ll5^ This result is remarkable 
when one recalls that there is no explicit bending rigidity. 

A plausible explanation 58] for the loss of the crumpled phase is that next-to-nearest neighbor excluded volume 
effects induce a positive bending rigidity, driving the model to the FLFP. The structure function of the self-avoiding 
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model has been computed numerically 59] and found to compare well with the analytical structure function for the 
flat phase of phantom fixed-connectivity membranes. In particular the roughness exponents are comparable. 

The induced bending rigidity may be lowered by taking a smaller excluded volume. |6flj The flat phase persists to 
very small values of a with eventual signs of a crumpled phase, probably due to effective loss of self-avoidance. A 
more comprehensive study, 61] in which the hard sphere radius is taken to zero with an excluded volume potential 
which is a function of the internal distance along the lattice, concluded that self-avoidance implies flatness in the 
thermodynamic limit of large membranes. 

Self-avoidance may also be implemented by modelling impenetrable triangular meshes. This has the advantage that 
there is no restriction on the bending angle between adjacent cells and therefore no induced bending rigidity. |62| 

The first simulations of the plaquette modell63| found a Hausdorff dimension in rough agreement in agreement with 
the Flory estimate 2.5 but this has not held up in subsequent work. A subsequent simulation 64J found dn ~ 2.3 and 
extensive recent work employing more sophisticated algorithms and extending to much larger membranes confirm the 
loss of the crumpled phase. (6^] 

Some insight into the lack of a crumpled phase for self-avoiding fixed-connectivity membranes is offered by the study 
of folding. [65l| ~ 171*1 Folding corresponds to the limit of infinite elastic constants 46] with the further approximation 
that the space of bending angles is discretized. One quickly discovers that the reflection symmetries of the allowed 
folding vertices forbid local folding (crumpling) of surfaces. There is therefore essentially no entropy for crumpling. 
There is, however, local unfolding and the resulting statistical mechanical models are non-trivial. The lack of local 
folding is the discrete equivalent of the long-range curvature-curvature interactions that stabilize the flat phase. The 
dual effect of the integrity of the surface (time-independent connectivity) and self-avoidance is so powerful that 
crumpling seems to be impossible in low embedding dimensions. 

2. The properties of the self- avoiding fixed point 




FIG. 13: The conjectured phase diagram for self-avoiding fixed-connectivity membranes in 3 dimensions. With any degree of 
self-avoidance the renormalization group flows are to the flat phase fixed point of the phantom model (FL) . 

For the physically relevant case d — i numerical simulations thus find that there is no crumpled phase. Furthermore, 
the flat phase is identical to the flat phase of the phantom membrane. '625 The roughness exponent C^sA from mmrerical 
simulations of self- avoidance at d = 3 using ball-and-spring models 72j and impenetrable plaquette models, 6 2j and 
the roughness exponent at the FLFF, Ea. H24() . compare extremely well 

CsA = 0.64(4) , C = 0.64(2) . (31) 



The numerical evidence thus strongly indicates that the SAFP is exactly the same as the FLFF and that the crumpled 
self-avoiding phase is absent in the presence of purely repulsive potentials (see Fig ll3|l . This conjecture is strengthened 
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by the finding that the Poisson ratio of self-avoiding membranes is the same as that of flat phantom membranes. ySj 
(see Sect II V|) . This identification of fixed points enhances the significance of the FLFP treated earlier. 

IV. POISSON RATIO AND AUXETICS 

In the classical theory of elasticity's^ an arbitrary deformation of a Z3-dimensional elastic body may be decomposed 
into a pure shear and a pure compression: 



l(Tr^)5,,] + -l 



{Tru)S.,,] + -{Tvu)6,, , (32) 



where Tr denotes the trace and the term in square brackets is a pure shear (volume-preserving but shape changing) 
while the second term is a pure compression (shape-preserving but volume-changing) . The elastic free energy is then 
given by 



Ui] - ^{Ti- u)Sij 



1 2 

+ -KiTTuf , (33) 



Fel = fJ- 

where fi is the shear modulus and K is the bulk modulus. This free energy may be written equivalently as 

1 2 

Fel = UijU^j + -A (Tr u) , (34) 
with the elastic Lame coefficient A related to the bulk and shear moduli by 

K^X+^. (35) 

For the physical membrane, D = 2, this reads K = X + fi. Thermodynamic stability requires that both K and fj, be 
positive, otherwise the free energy could be spontaneously lowered by pure compressional or pure shear deformations, 
respectively. 

The Poisson ratio a is defined as the ratio of transverse contractile strain to longitudinal tensile strain for an elastic 
body subject to a uniform applied tension T. For tension applied uniformly in, say, the j:-direction 

a = -j—j- > (36) 
ox/x 



the Poisson ratio is easily found to be 



K + n 



- - \ ) (38) 



for D = 2, and 



2 V 3/v + M 

for D = 3. Thermodynamic stability is only possible for —1 < cr < 1 for D = 2 and —1 < cr < iforZ? = 3. The upper 
bounds (1 and 1/2 respectively) are approached for materials that have vanishing shear modulus compared to their 
bulk modulus (rubber-like) and the lower bounds (-1) for materials with negligible bulk modulus in comparison to their 
shear modulus ( "anti-rubber" ) [t^ . Clearly, the Poisson ratio may be negative (auxetic) for if < ^ (D=2) and K < ^fi 
(D=3). Most materials get thinner when stretched and fatter when squashed - auxetic materials are uncommon. The 
earliest known example, dating from more than a century ago, is that of a pyrite (FeS2) crystal 75], which has a 
Poisson ratio, in certain crystallographic directions, of cr « —0.14. More recently, some isotropic polyester foams have 
been created with Poisson ratios as large as cr '--^ — O.tItJIt^. The potential of auxetic materials in materials science 
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is nicely reviewed in Ref.[79j. One of the rare naturally occurring auxetics is Si02 in its a-crystobalite phase. |8Cl l8l| 
Cristobalite is one of the three distinct crystalline forms of Si02, together with quartz and tridymite. Its Poisson ratio 
reaches a maximum negative value of —0.5 in some directions, with orientationally-averaged values for single-phased 
aggregates of —0.16. 

The underlying mechanism driving fixed-connectivity membranes auxetic (cr = —1/3) has schematic similarities 
to that illustrated in Fig ll4l Submitting a membrane to tension will suppress its out-of- plane fluctuations, forcing 
it entropically to expand in both in-plane directions. More physically, the out-of-plane undulations renormalize the 
elastic constants (the Lame coefficients), in such a way that the long- wavelength bulk modulus is less than the 
shear modulus, which is the signature of a two-dimensional auxetic material. The soft matter origin of the universal 
negative Poisson ratio of fixed-connectivity membranes provides a fundamentally new paradigm for the design of novel 
materials. The best current experimental measurements of the Poisson ratio of the red blood cell cvtoskeleton [s^ 
find a w +1/3 from separate determinations of the bulk and shear modulus. The cytoskeleton still has the fluid 
lipid bilayer attached, however, and this may influence the pure cytoskeletal elasticity. A direct measurement of the 
Poisson ratio for a flexible flxed-connectivity membrane remains an important and challenging task. 




FIG. 14: Mechanical model of an auxetic material: (a) in the absence of applied stress and (b) under applied lateral stress T . 
The lateral stretching accompanying the applied stress forces the material out in the transverse dimension. 



Auxetic materials have desirable mechanical properties such as higher in-plane indentation resistance, transverse 
shear modulus and bending stiffness. They have clear applications as sealants, gaskets and fasteners. They may also 
be promising materials for artificial arteries, since they can expand to accommodate sudden increases in blood flow. 

We can model a realistic fixed-connectivity membrane with an elastic free energy and either large bending rigidity or 
self- avoidance. This is of practical importance in modelling since, for example, we may replace the more complicated 
non-local self-avoidance term with a large bending rigidity. 

It would be very interesting to know if nature utilizes the auxetic character of the red-blood cell spectrin cytoskeleton 
in the elastic deformations of red blood cells as they pass through fine blood capillaries. As such cells deform, the 
membrane skeleton can unfold, which might help to transport large molecules or expose reactive chemical groups, [s^ 



V. ANISOTROPIC MEMBRANES 



An anisotropic membrane is a fixed-connectivity membrane in which the elastic moduli or the bending rigidity 
in one distinguished direction are different from those in the remaining D — \ directions. Such a membrane may 
be described by a d-dimensional vector r(x^,j/), where now the D dimensional internal coordinates are split into 

— 1 XX coordinates and the orthogonal distinguished direction y. 

Requiring invariance under translations, 0(d) rotations in the embedding space and 0(D — 1) rotations in the 
internal space, the equivalent of Eq.Q becomes 
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+K^ydlr- dlf+ t^id^ry + tyidyry 
+ ^{d;tr- d^rr + ^{dyr- dyv)' 

+ ^ J d^x J d^x'(5'^(f(x) - f(x')) . (39) 

This model has eleven free parameters - three distinct bending rigidities, K±,Ky and Kj_y, seven elastic moduli, 
t±,ty, Uyy, u ±y , v±± Slid V ±y , Slid tlic strcHgth of sclf" avoidaucc coupling b. 
As before we decompose displacements as 



f(x) = (c±x_L + u^(x), CylJ + Uyix.), K(x)j , (40) 

with being the D — 1-dimensional intrinsic phonon modes, Uy the intrinsic phonon mode in the distinguished 
direction y and h the d — D-dimensional out-of-plane fluctuation mode. If C-L = Cy = 0? the membrane is crumpled 
and if both (± and Cy do not vanish the membrane is flat. There is, however, the possibility that C± = ^^^d 
or 7^ and Cy = 0. This describes a tubular phase, in which the membrane is crumpled in some internal directions 
but flat in the remaining ones.[8^ l85j. Fig ll5l displavs a typical equilibrium conflguration from the tubular phase, 
along with the low and high-temperature flat and crumpled phases for a phantom anisotropic membrane. 

Let's deal with the phantom anisotropic membrane flrst. Both analvtical|8fi|| and numerical work|^ has established 
that the phase diagram contains a crumpled, tubular and flat phase. The crumpled and flat phases are equivalent 
to the isotropic ones, so anisotropy turns out to be an irrelevant interaction in those phases. The new physics is 
contained in the tubular phase. 



A. Phantom Tubular Phase 



1. The Phase diagram 

We first describe the mean field theory phase diagram and then the effect of fluctuations. There are two situations 
dep ending on the value of a certain function A, which depends on the elastic constants u±±,v±y,Uyy and v i i .|84L 

For A > the mean field solution exhibits crumpled, flat and tubular phases. When > and t_L > the model 
is crumpled. Lowering the temperature so that one of the t couplings becomes negative drives the membrane to the 
tubular phase (either a _L or y-tubule). Lowering the temperature still further flattens the membrane. For A < 
the flat phase disappears from the mean field solution, leaving only the crumpled and tubular phases separated by a 
continuous transition. Tubular phases are the stable low temperature stable phases in this regime. This mean field 
result is summarized in Fig llBI 

Beyond mean field theory, the Ginzburg criterion applied to this particular model suggests that the phase diagram is 
stable for physical membranes Z? = 2 at any embedding dimension d. The mean field description should be qualitative 
correct even for the full model. 

Numerical simulations have spectacularly confirmed this beautiful analytic prediction, [s^ Changing the temperature 
generates a sequence of continuous phase transitions crumpled-to-tubular and tubular-to-flat, in total agreement with 
the A > case above (see Fig [TC|) . 

2. The Crumpled Anisotropic Phase 
In this phase ty > and t± > 0, and the free energy Ea. H39() reduces , for D >2, to 

F{f{x)) = ^ / d^-^x±dy [t±{d^r)^ + ty{dyff] + Irrelevant Terms . (41) 
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The three phases of anisotropic 
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FIG. 16: The phase diagram for anisotropic phantom membranes. 



By redefining the y coordinate to be y' — j-y this reduces to Eq.®, with t = t±. Anisotropy is clearly irrelevant in 
the crumpled phase. 



3. The Flat Phase 



In the flat phase intrinsic anisotropics are only apparent at short-distances and therefore should be irrelevant in 
the infrared limit. This argument may be made more precise. [s^ The flat phase is thus equivalent to the flat phase 
of isotropic membranes. 
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B. The Tubular Phase 

We now turn to the study of the novel tubular phase, both in the phantom case and with self-avoidance. Since the 
physically relevant case for membranes is D — 2 the y-tubular and _L-tubular phase are the same. 

The key critical exponents characterizing the tubular phase are the size (or Flory) exponent i/, giving the scaling 
of the tubular diameter Rg with the extended (Ly) and transverse {L±) sizes of the membrane, and the roughness 
exponent ( associated with the growth of height fluctuations hrms (see Fig ll7f) : 



Rg{L_L,Ly) CX L'^SB.{Ly/L]_) ; 

hrms{L±, Ly) cx L^ySh{Ly / L^^_) , 

where Sb. and Sh are scaling functions |84L l85l| and z is the anisotropy exponent. 



(42) 




FIG. 17: A schematic illustration of a tubular configuration indicating the radius of gyration Rg and the height fluctuations 

hrms • 

The general free energy described in Ea. lp!(^ may be simphfied considerably in a y-tubular phase. [s^ l9(Tj|: 



F{u,h) = 



-g_i_{daU + dahdyh) 
1 



+ 9yidy^+^idyhfr 



(43) 



reducing the number of free couplings to five. The coupling g±, furthermore, is irrelevant by standard power counting. 
The most natural assumption is to set it to zero. In that case the phase diagram one obtains is shown in Fig ll8l 
Without self-avoidance, i.e. b = 0, the Gaussian Fixed Point (GFP) is unstable and the long-wavelength behavior of 
the membrane is controlled by the tubular phase fixed point (TPFP). Any amount of self-avoidance, however, leads 
to a new fixed point, the Self-avoiding Tubular fixed point (SAFP), which describes the large distance properties of 
self-avoiding tubules. 

Radzihovsky and Toner advocate a different scenario. |85j For sufficiently small embedding dimensions d, including 
the physical d = 3 case, these authors suggest the existence of a new bending rigidity renormalized fixed point (BRFP), 
which is the infra-red fixed point describing the actual properties of self-avoiding tubules (see Fig. [T^ . 

Here we follow the arguments presented in Refs.|89l.l90l| and consider the model defined by Ea. H43|) with the g^-term 
set to zero. One can prove then than there are some general scaling relations among the critical exponents. All three 
exponents may be expressed in terms of a single exponent 



3 1-D 



2 2z 
V — C^z . 



(44) 



Remarkably, the phantom case, as described by Ea. (|43|l . can be solved exactly. The result for the size exponent is 



5 -2D 1 
Vvh(P) = — ^ — , i'pft(2) = - 



(45) 
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FIG. 18: The phase diagram for self-avoiding anisotropic membranes with the Gaussian fixed point (GFP), the tubular phase 
fixed point (TPFP) and the self- avoidance fixed point (SAFF). 




FIG. 19: The phase diagram for self-avoiding anisotropic membranes with the Gaussian fixed point (GFP), the tubular phase 
fixed point (TPFP), the self- avoidance fixed point (SAFP) and the bending rigidity fixed point (BRFP). 

with the remaining exponents following from the scaling relations Eg . (|44|) . 

The self-avoiding case may be treated with techniques similar to those in isotropic case. The size exponent may be 
estimated within the Flory approximation, yielding 



UFi = 2/dH = -—■ ■ (46) 
a+1 

The Flory estimate is an uncontrolled approximation. Fortunately, a e-expansion, adapting the multi-local operator 
product expansion technique 91]^ 93J to the case of tubules, is also possible. 89, 90] The resulting rcnormalization 
group /3-functions provide evidence for the phase diagram shown in Fig ll8l Extrapolation techniques also provide 
estimates for the size exponent, the most accurate value being v ~ 0.62 for the physical case. The rest of the exponents 
may be computed from the scaling relations. 

Numerical simulations so far, however, do not find a tubular phase in the case of strict self- avoidance. 

VI. ORDER ON CURVED SURFACES 

Imagine we instantaneously freeze a fluctuating membrane so that it has some fixed but curved shape. We can then 
ask about the nature of the ground state of particles distributed on this surface and interacting with some microscopic 
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pair-wise repulsive potential. The relevant physics is clearly related to the infinite bending rigidity limit (flat phase) 
of elastic membranes. In such a membrane the topology and topography are fixed. 




FIG. 20: A 2.5 micron scan of 0.269 micron diamet er polystyrene spheres crystallized into a regular triangular lattice 
from |http://invsee.asu.edu/nmodules/spheresmod/| 



taken 



Spherical particles on a fiat surface pack most efficiently in a simple triangular lattice, as illustrated in Fig l2UI In the 
dense limit each particle "kisses" six of its neighbors. '94] Such six-coordinated triangular lattices cannot, however, be 
perfectly wrapped on the curved surface of a sphere; topology alone requires there be defects in coordination number. 
The panels on a soccer ball and the spherical carbon molecule Cgo (buckyball) [qtL l9^ are good illustrations of the 
necessity of defects for a spherical triangulation ~ they have 12 pentagonal faces (each the dual of a 5-coordinated 
defect) in addition to 20 hexagonal faces (each dual to a regular 6-coordinated node.) The necessary packing defects 
can be characterized by their topological charge, q, which is the departure of their coordination number c from the 
preferred fiat space value of 6 (g = 6 — c). These coordination number defects are point-like topological defects 
called d isclin ations 99] and they detect intrinsic Gaussian curvature located at the defect. A profound theorem 
of Euler p'ool IToi| states that the total disclination charge of any triangulation whatsoever of the sphere must be 
12! 1q2 a total disclination charge of 12 can be achieved in many ways, however, which makes the determination 
of the minimum energy configuration of repulsive particles, essential for crystallography on a sphere, an extremely 
difficult problem. This was recognized nearly 100 years ago by J.J. Thomson [l03j . who attempted, unsuccessfully, to 
explain the periodic table in terms of rigid electron s hells. Similar problems arise in fields a s diverse as mult i-electron 
bubbles in supcrfluid hefiu m. llQ^I virus morphology. [lOfij ~ |l08l | protein s-lavers. [l09l IllOj l giant molecules |l 111 HT^ 
and information processing. |95ll96l | Indeed, both the classic Thomson problem, which deals with particles interacting 
through the Coulomb potenti al, a nd its generalization to other interaction potentials, are still open problems after 
almost 100 years of attention. |ll3j ~ |115| 

The spatial curvature encountered in curved geometries adds a fundamentally new ingredient to crystallography 
not found in the study of order in spatially flat systems. As the number of particles on the sphere grows, isolated 
charge 1 defects (5s) will induce too much stra in. This s train can be relieved by introducing additional dislocations, 
consisting of pairs of tightly bound 5-7 defects flT^ IllTj . which don't spoil the topological constraints because their 
net disclination charge is zero. Dislocations, which are themselves point-like topological defects in two dimensions, 
disrupt the translational order of the crystalline phase but are less disruptive of orientational order. [TT^ 

Recent work on an experimental realization of the generalized Thomson problem has allowed us to explore the lowest 
energy configuration of the dense pack ing of repulsive particles on a spherical surface and to confront a previously 
developed theory with experiment • |l 1 8| . We create two-dimensional packings of colloidal particles on the surface of 
spherical water droplets and view the structures with optical microscopy. Above a critical system size, the thermally 
equilibrated colloidal crystals display distinctive high-angle grain boundaries, which we call "scars" . These grain 
boundaries are found to end entirely within the crystal, which is never observed on fiat surfaces because the energy 
penalty is too high. 

The experimental system is based on the self-assembly of one micron diameter cross-linked polystyrene beads 
adsorbed on the surface of spherical water droplets (of radius R), themselves suspended in a density-matched oil 
mixture. |ll9l | The polystyrene beads are almost equally happy to be in oil or water (the bead/oil surface tension is 
close to the bead/water surface tension) and therefore diffuse freely until they find the oil- water interface and stick 
there. Particle assembly on the interf ace o f two distinct liquids dates to the pioneering work of Pickering |l20j and 
was beautifully exploited by Pieranski |l2l) some time ago. The particles are imaged with phase contrast using an 
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FIG. 21: Light microscope images of particle-coated droplets. Two droplets (A) and (C) are shown, together with their 
associated defect structures (B) and (D). Panel (A) shows an ~ 13% portion of a small spherical droplet with radius R — 12.0 
microns and mean particle spacing a = 2.9 microns (R/a = 4.2), along with the associated triangulation (B). Charge +1( — 1) 
disclinations are shown in red and yellow respectively. Only one +1 disclination is seen. Panel (C) shows a cap of spherical 
colloidal crystal on a water droplet of radius R — 43.9 microns with mean particle spacing a = 3.1 microns (R/a = 14.3), along 
with the associated triangulation (D). In this case the imaged crystal covers about 17% of the surface area of the sphere. The 
scale bars in (A) and (C) are 5 microns. 



We analyze the lattice configurations of a collection of 40 droplets. A typical small spherical droplet with system 
size, R/a = 4.2, where a is the mean particle spacing, is shown in Fie:.l21lA.. The associated Delaunay triangulation is 
shown in Fig. l2Hi. The only defect is one isolated charge +1 disclination. Extrapolation to the entire surface of the 
sphere is statistically consistent with the required 12 total disclinations. 

Qualitatively different results are observed for larger droplet sizes as defect configurations with excess dislocations 
appear. Although some of these excess dislocations are isolated, most occur in the form of distinctive (5 — 7— 5 — 7 — 
• • • — 5) chains, each of net charge +1, as shown in Fig. l21L ). These chains form high-angle (30°) grain boundaries, or 
scars, which terminate freely within the crystal. Such a feature is energetically prohibitive in equilibrium crystals in 
flat space. Thus, although grain boundaries are a common feature of 2D and 3D crystalline materials, arising from a 
mismatch of crystallographic orientations across a boundary, they usually terminate at the boundary of the sample in 
flat space because of the excessive strain energy associated with isolated terminal disclinations. Termination within 
the crystal is a feature unique to curved space. 

Of key interest is the number of excess dislocations per chain as a function of the dimensionless system size R/a. 
This is plotted in Fig. [53 Scars only appear for droplets with R/a > 5. These results provide a critical confirmation 
of a theoretical prediction that R/a must exceed a threshold value (R/a)e ~ 5 , corresponding to M ~ 360 particles, 
for excess defects to proliferate in the ground state of a spherical crvstal. [l23j The precise value of (R/a)c depends 
on details of the microscopic potential, but its origin is easily understood by considering just one of the 12 charge -1-1 
disclinations required by the topology of the sphere. In fl at space such a topological defect has an associated energy 
that grows quadratically with the size o f the system. [ll7j since it is created by excising a 27r/6 wedge of material and 
gluing the boundaries together. [ill Ill7j The elastic strain energy associated with this defect grows as the area. In the 
case of the sphere the radius plays the role of the system size. As the radius increases, isolated disclinations become 
much more energetically costly. This elastic strain energy may be reduced by the formation of linear dislocation 
arrays, i.e. grain boundaries. The energy needed to create these additional dislocation arrays is proportional to a 
dislocation core energy Ec and scales linearly with the system size. 123] Such screening is inevitable in flat space (the 
plane) if one forces an extra disclination into the d e fect- free ground state. Unlike the situation in flat space, grain 
boundaries on the sphere can freely terminate. |l23j ~ |126| consistent with the scars seen on colloidal droplets. 

A powerful a nalytic approa ch to determining the ground state of particles distributed on a curved surface has 
been developed. |l23l Il26l Il27j The original particle problem is mapped to a system of interacting disclination defects 
in a continuum elastic curved background. The defect-defect interaction is universal with the particle microscopic 
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FIG. 22: Excess dislocations as a function of system size. The number of excess dislocations per minimal disclination A'^ as a 
function of system size R/a, with the linear prediction given by theory shown as a solid red line. 



potential determining two free parameters - the Young modulus Kq of the elastic background and the core energy Ec 
of an e leme ntary disclination. A rigorous geometrical derivation of the effective free energy for the defects is given 
in Ref. |l28l | An equivalent derivation may also be given by integrating out the phonon degrees of freedom from the 
elastic Hamiltonian,[23 with the appropriate modifications for a general distribution of defects. The energy of a 
two-dimensional crystal embedded in an arbitrary frozen geometry described by a metric gij (x) is given by 

H = Eo +— / / da{x)d(j{y) 



2 _ 

[s{x)-K{x)]^[{s{y)-K{y)] \ , (47) 



where the integration is over a fixed surface with area element da^x) and metric g.ij, K is the Gaussian 



curvature, Y is the Young modulus in flat space and s{x) = 'Y^f^i'^cii5{x,Xi) is the disclination density 



6{x,Xi) — 5{x — Xi) / det{gij) . Here 5- and 7-fold defects correspond to qi = +1 and —1, respectively. Defects 

like dislocations or grain boundaries can be built from these N elementary disclinations. i?o is the energy correspond- 
ing to a perfect defect-free cryst al wi th no Gaussian curvature; Eq would b e th e ground state energy for a 2D Wigner 
crystal of electrons in the plane.|l29l| Ea. l|47|l . restricted to a sphere, gives |l2,'ij 

-Y w 

i=l j=l 

+ NE, , (48) 

where Ec is a defect core energy, R is the radius of the sphere and x is a function of the geodesic distance Pij between 
defects with polar coordinates {6^,tp'^ ; 9^,tp^): 

In z 



XW = 1 + / dz . (49) 

Jo 

The potential is attractive for opposite charged defects and repulsive for like-charged defects. Many predictions of this 
model are universal in the sense that they are insensitive to the microscopic potential. This enables us to make definite 
predictions even though the colloidal potential is not precisely known. It also means that our model system serves as a 
prototype for any analogous system with repulsive interactions and spherical geometry. To further test the validity of 
this approach, we show a typical ground state for large M in Fig.|2Sl The system size here is R/a = 12, similar to the 
droplet in Fig.l21L). The results are remarkably similar to the experimentally observed configuration in Fig.l21U: the 
only difference is a result of thermal fluctuations, which break the two defect scars in the experiment. This agreement 
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FIG. 23: Model grain boundaries. This image is obtained from a numerical minimization for a system size comparable to the 
large droplet in Fie l21f c.d'l. 




FIG. 24: The polyhedral siliceous cytoskeleton of the unicellular ocean organism Aulosphaera. 



between theory and experiment also provides convincing evidence that these scars are essential components of the 
equilibrium crystal structure on a sphere. 

The theory predicts that an iso lated charge +1 disclination on a sphere is screened by a string of dislocations of 
length cos~^(5/6)i? ~ 0.59f?. [T23j One can use the variable linear density of dislocations to compute the total number 
of excess dislocations iV in a scar. One finds that TV grows for large (R/a) as f [711 - 5 cos-^S/e)] f w 0.4lf, 
independently of the microscopic potential. This prediction is universal, and is in remarkable agreement with the 
experiment, as shown by the solid line in Fig. 

We expect these scars to be widespread in nature. They should occur, and hence may be exploited, in sufficiently 
large stiff viral protein capsids, giant spherical fuUerenes, spherical bacterial surface layers (s-layers), provided that 
the spherical geometry is not too distorted. Terminating strings of he ptago ns and pentagons might serve as sites for 
chemical reactions or even initiation points for bacterial cell division[l09j and will surely influence the mechanical 
properties of spherical crystalline shells. 

The polyhedral silice ous cytoskeleton of the unicellular non-motile ocean organism Aulosphaera (a member of the 
species Phaeodaria \l3d^ ) is shown in Fig l24l A triangulation revealing three scars, two of which are branched, is shown 
in Fig[2Hl The skeleton itself is such a perfect triangular lattice that it coincides with the Delaunay triangulation 
determined by its vertices. The case of viral capsids has been analyzed in Ref. 108], where it is shown that, rather 
than scarring, icosahedral packings become unstable to faceting f or su fficiently large virus size, in analogy with the 
buckling instability of disclinations in two-dimensional crystals. Il3l| 
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FIG. 25: The Delaunay triangulation of the Aulosphaera cytoskeleon above. 




FIG. 26: A typical ground state for a toroidal hexatic. Five-fold disclinations are shown as solid circles (red) and 7-fold 
disclinations as solid squares (blue). 



Scarred spherical crystals may provide the building blocks (atoms) of micron-scale molecules |l32llT33| and materials. 
While topology dictates the overall number of scars (12), the details of the geometry and defect energetics determine 
the length and structure of the scars themselves. It is possible that scarred co Uoid osomes will ultimately yield complex 
self- assembled materials with novel mechanical or optoelectronic properties. 51325 

New structures arise if one changes the structure o f the colloid of the topology of the surface they coat. Nelson has 
analyzed the case of nematic colloids coating a sphere. |l33j | In this case the preferred number of elementary disclination 
defects is 4, allowing for the possibility of colloidal atoms with tetrahedral functionality and sp^-typ e bo nding. The 
case of toroidal templates with 6-fold bond-orientational (hexatic) order has been analyzed recently. |l27j It is found 
that defects are energetically favored in the ground state for fat torii or moderate vesicle size. A schematic of a 
"typical" ground state is shown in FigEHl 
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